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Abstract 

The naive time reversal odd ("T-odd") parton distribution hj~, the so-called Boer- 
Mulders function, for both up (u) and down (d) quarks is considered in the diquark 
spectator model. While the results of different articles in the literature suggest that the 
signs of the Boer-Mulders function in semi-inclusive DIS for both flavors u and d are the 
same and negative, a previous calculation in the diquark-spectator model found that w ' 
and hf^ have different signs. The flavor dependence is of significance for the analysis of 
the azimuthal cos(2</>) asymmetries in unpolarized SIDIS and DY-processes, as well as for 
the overall physical understanding of the distribution of transversely polarized quarks in 
unpolarized nucleons. We find substantial differences with previous work. In particular 
we obtain half and first moments of Boer-Mulders function that are negative over the full 
range in Bjorken x for both the u- and d- quarks. In conjunction with the Collins function 
we then predict the cos(20) azimuthal asymmetry for ir + and ir~ in this framework. We 
also find that the Sivers u- and d- quark are negative and positive respectively. As a by- 
product of the formalism, we calculate the chiral-odd but "T-even" function h^ L , which 
allows us to present a prediction for the single spin asymmetry A s ™^ 2 ^ for a longitudinally 
polarized target in SIDIS. 



1 Introduction 

Naive time reversal-odd ("T-odd") transverse momentum dependent (TMD) parton distribu- 
tions (PDFs) have gained considerable attention in recent years. Theoretically it is expected 
that they can account for non-trivial transverse spin and momentum correlations such as single 
spin asymmetries (SSA) in hard scattering processes when transverse momentum scales are 
on the order of intrinsic transverse momentum of quarks in hadron, namely Pt ~ k± -C \J Q 2 . 
Experiments are being performed [H El E] and proposed [U [5] to test these hypotheses by 
measuring transverse SSAs (TSSAs) and azimuthal asymmetries (AAs) in hard scattering 
processes such as semi-inclusive DIS (SIDIS) or the Drell-Yan process (DY). A prominent 
example of such a "T-odd" PDF is the Sivers function f^f (6j [7] which explains the observed 
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SSA in SIDIS for a transversely polarized proton target by the HERMES collaboration [JJ. 
It correlates the intrinsic quark transverse momentum and the transverse nucleon spin. The 
corresponding SSA on a deuteron target measured by COMPASS [2] vanishes, indicating a 
flavor dependence of the Sivers function. 

Another leading twist "T-odd" parton distribution, introduced in [8], correlates the trans- 
verse spin of a quark with its transverse momentum within the nucleon, the so called the 
Boer-Mulders function h^. It describes the distribution of transversely polarized quarks in an 
unpolarized nucleon. 

Theoretically, twist two "T-odd" PDFs are of particular interest as they formally emerge 
from the gauge link structure of the color gauge invariant definition of the quark-gluon-quark 
correlation function [9j QUI [UJ. This gauge link not only ensures a color gauge invariant 
definition of correlation functions, but it also describes final state interactions [32] and initial 
state interactions [13] which are necessary to generate SSA [12j [14j [TT] . Assuming factorization 
of leading twist SIDIS spin observables in terms of the "T-even"pI] and "T-odd" TMD PDFs 
and fragmentation functions (FFs), Ref. [8] shows how spin observables in SIDIS can be 
expressed in terms of convolutions of these functions. Formal proofs of factorization of leading 
twist SIDIS spin observables were presented later in Refs. [T6l [TTl fT8] . 

Apart from the leading twist transverse SSA, measurements were also performed in SIDIS 
on sub-leading SSA (i.e. they are suppressed like 1/Q, where Q is the virtuality of the ex- 
changed photon). In particular, the asymmetry for a longitudinally polarized target was 
measured by HERMES [IH [20l EU, |22] , whereas a non- vanishing beam-spin asymmetry was 
reported by CLAS [23j[24]. It was shown that final state interactions contribute also to these 
types of single-spin asymmetries [25j [261 EH EE] • Subsequently, this effect was described by 
the introduction of heretofore unknown sub-leading twist "T-odd" PDFs [29J; a complete list 
of these PDFs was presented in Ref. [30]. These sub-leading twist "T-odd" PDFs discovered 
in this work were then incorporated into the tree-level formalism |31] completing the original 
work of [15]. 

In this paper we focus on the flavor dependence of the leading twist-2 "T-odd" parton 
distributions in semi-inclusive DIS, i.e. Boer-Mulders function hj~, which is also chirally odd, 
and the Sivers function f^p (keeping in mind that "T-odd" parton distributions in the Drell- 
Yan process flip their sign [9]). The Boer-Mulders function is particularly important for the 
analysis of the azimuthal cos(20) asymmetry in unpolarized SIDIS and Drell-Yan. While in a 
partonic picture of the unpolarized cos(20) asymmetry in SIDIS, the Boer-Mulders function 
is convoluted with the "T-odd" (and chiral-odd) Collins fragmentation function [32], the 
corresponding cos(2</>) asymmetry in DY includes a convolution of the type hj^'Sihi [33] (where 
hi is the Boer-Mulders function for anti-quarks). Although these azimuthal asymmetries 
were measured in SIDIS by the ZEUS collaboration [34l [35] and in DY [361 [371 [38], little 
is known about the Boer-Mulders function. Of particular interest is the sign for different 
flavors u and d since this significantly affects predictions for these asymmetries (see Ref. 
[39]). The flavor dependence of h± was studied in the MIT-bag model [40] as well as in a 
spectator diquark model [UJ, and a large N c analysis of TMDs was performed in Ref. [42] . 
Model calculations of chirally odd generalized parton distributions (GPDs) [43] and a study 
of generalized form factors in lattice QCD [44] give indications about the flavor dependence 
of by means of non-rigorous and model-dependent relations between GPDs and transverse 
momentum dependent PDFs which were proposed and discussed in Refs. [45[ [46], [47] . All of 
these theoretical and phenomenological treatments suggest an equal (and negative) sign for 
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the Boer-Mulders function for both u- and d-quarks, with the exception of the calculation 
in the diquark spectator model which results in opposite signs for u and d. The purpose of 
this paper is to consider the flavor dependence of h± and extend our earlier work on this 
subject [48j[49]. Additionally, we consider the flavor dependence of the "T-even" function, h^ L 
which is also of interest in exploring the transverse momentum and quark spin correlations in 
a longitudinally polarized target |50J. 



2 "T-Odd" PDFs in the Spectator Model 

Transverse momentum quark distribution and fragmentation functions contain essential non- 
perturbative information about the partonic structure of hadrons. Practically speaking their 
moments are calculable from first principles in lattice QCD. A great deal of understanding 
has also been gained from model calculations using the spectator framework. In addition 
to exploring the kinematics and pole structure of the TMDs [5TJ [52j [18] phenomenological 
estimates for parton distributions [53] and fragmentation functions p3] for "T-even" PDFs 
and for "T-odd" PDFs [121 H31 SH UHl HH [55] EH HZj have been performed. We extend 
these studies to explore the flavor dependence of the "T-odd" pdfs adopting the factorized 
approach used in Refs. [541 1531 149] , 

We start (cf. [53]) from the definition of the fully- unintegrated, color gauge invariant, 
quark-quark correlator 

Syfo P,S) = Y,J 704 e w< (P, S\ ^(0) W[0 | oo, 0, T ] \X)(X\ W[oo, | £] \P, S), 

(1) 

where the gauge link indicated by the (straight) Wilson line is given by 

b 

(2) 



W[a\b] = ^exp|-i 5 ^ WA^s) 



In an arbitrary gauge there is a Wilson line at light cone infinity pointing in transverse direc- 
tions [TOldl]. Here, we work in Feynman gauge where the transverse Wilson line vanishes [10J. 
In the definition ([!]) we insert a complete set of intermediate states 1 = ^2 X \X)(X\. In the 
diquark model the sum over a complete set of intermediate on-shell states \X) is represented 
by a single one-particle diquark state | dq; pd q , A), where pd q is the diquark momentum and A 
its polarization. Since the diquark is "built" from two valence quark it can be a spin particle 
(scalar diquark) or a spin 1 particle (axial- vector diquark). By applying a translation on the 
second matrix element in Eq. ((TJ) we can integrate out £, perform the momentum integration 
over the diquark momentum pd q , and obtain 

*ijfaP,S) = 6{{P ~ ~ m ; ] 3 &{P0 - P ° ] (P, S\ & (0) W[0 | oo, 0, T ] \dq-P- p, A) x 

(dq; P-p,A|W[oo,0,0V|0]^(0)|P,5>. (3) 

The essence of the diquark spectator model is to calculate the matrix elements in Eq. J3]) by 
the introduction of effective nucleon-diquark-quark vertices. 
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Figure 1: Different vertices for the axial- vector diquark. Left Panel: Nucleon-diquark-quark vertex. 
Right Panel: Diquark-gluon vertex. 



For "T-even" parton distributions such as the unpolarized PDF f\ one obtains a non- 
vanishing result at leading order (in the nucleon-diquark-quark coupling) with a trivial con- 
tribution from the Wilson line, i.e. at tree level. In this case the matrix element (dq\ip\P) is 
depicted in the Left Panel of Fig. CD For a scalar and an axial-vector diquark different ver- 
tices have to be chosen. The most general nucleon-diquark-quark vertices for off-shell particles 
were presented in Ref. [56]. For the matrix elements (dq\ip\P), the nucleon is on-shell which 
reduces the amount of structures of the vertices of Ref. [56J. In the following we work with 
the nucleon-diquark-quark vertices which were used in Ref. [53j to compute "T-even" PDFs. 
They read for a scalar and an axial-vector diquark 



T s (N)=g sc (p 2 ) ; r» x (N) 



9ax(p 2 

V3 



"75 



9 M 



(4) 



g(p 2 ) are form factors depending on the quark momentum p. They are introduced to yield a 
more realistic description of the non-perturbative nature of the quark-quark correlator, and 
are specified below. R g is a ratio of coupling constants, since both structures in the nucleon- 
( axial- vector) diquark-quark coupling can in principle have different couplings. 

To leading order, the matrix elements are given by the following expressions for a scalar 
and axial vector diquark, 



(sdq; P-p\il>i(p)\P,S) 



{adq;P-p;\\^(0)\P,S) 



ig S c(p 



2 ^ + m q )u(P,S)] t 
p 2 — m 2 + iO 



m, 



•)-*> [r 



R 9lP\ U 



p 2 — m 2 + iO 



(5) 



(6) 



where the polarization vector of the axial-vector diquark is given by e^, u(P,S) denotes the 
nucleon spinor, and M and m q are nucleon and quark masses, respectively. In this paper we 
consider the diquark as a particle with mass m s , and the polarization sum for the axial- vector 
diquark is 



J2 £ U P ~ p ; A ) £ ^ P - A ) = -9v» + 



(P-p)^P-p\ 



mi 



(7) 



The unpolarized TMD /i is obtained by inserting Eqs. JH) and ((H) into Eq. J3|) and projecting 
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fi from the quark-quark correlator (see e.g. [HHl EI 



2fi(x,f T ) = \j 'dp' (Tr [ 7 +$(p;P,5)] +Tr [ 7 +$(p; P, -S)]) 



(8) 



where the "+" sign of the 7-matrix denotes the usual light cone component (a^ = l/v2(a° db 
a 3 )). The results for /1 in the scalar and axial vector diquark sectors read 

/f(^t) = tJ~^ \9sc(p 2 ) | 2 r ^ " ^ 2 + (xM + m g ) 2 l , (9) 



m 2| 



(10) 



where rh 2 = xm 2 — x(l — x)M 2 + (1 — x)m 2 . To shorten the notation we introduce a function 
lZf x depending on x and pr, and the model parameters R g and the set of masses, i.e. {Ai} = 
{M, m s , m q }, to be fixed below. 

Another "T-even" function of interest is the distribution of transversely polarized quarks 
in a longitudinally polarized target, 

2\ P pi r hi L (x,f T ) = y J *P~ ( Tr [tVts^Cp; P, S L )] - Tr [<y + <y i <y B *(l>\ P, Sl)] } , (11) 
where Xp is the target helicity, and Sl is the spin 4- vector in longitudinal direction, i.e. 



Sl 



Ap p- Xp p+ q 
M r ' M ' UT 



. By applying the same methods as for /1, we obtain 

sc (p 2 )\ 2 (l-x)M(xM 
( 27r ) 3 [p% + m 2 ] 



,±, sc , -ax \g S c(p 2 )\ 2 (l-x)M(xM + m q ) 

n lL \ x >Pt) — 77T^3 7^> T7T2 ' \ 1Z ) 



*iH*.*> - ^^^U.,^ • (13) 

where for brevity 7£" x and Tlf L ax are given in Appendix O 

By contrast, "T-odd" PDFs cannot be generated by simply considering the tree-level di- 
agram in the Left Panel of Fig. [TJ In the spectator framework the "T-odd" PDFs [12] are 
generated by the gauge link in Eq. ([T]) [HI 33 HH [TS] . Thus, the leading contribution can be 
obtained by expanding the exponential of the gauge link up to first order. This contribution 
results in a box diagram as shown in the Left Panel of Fig. [21 which contains an imaginary 
part necessary for "T-odds". We restrict ourselves to the case where one gluon models the final 
state interactions. The contribution of the gauge link is represented in the Left Panel of Fig. [2] 
by the double (eikonal) line and the eikonal vertex yielding a contribution to the box diagram 

x {-ie q v x ), (14) 



[l-v + iO] 



where I is the loop momentum, e q the charge of the quark and v is a light cone vector represent- 
ing the direction of the Wilson line. In order to evaluate the box diagram we need to specify 
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the gluon-diquark coupling. With an one-gluon exchange approximation in mind we use the 
gluon-diquark coupling for a scalar diquark, and for an axial-vector diquark we use a general 
axial-vector-vector that models the composite nature of the diquark through an anomalous 
magnetic moment k [57]. In the notations of Fig. [1] (Right Panel) the gluon-diquark vertices 
read 

r» = -ie dq ( Pl + P2 r, (15) 

I™ = -ie dq [g^( Pl +p 2 r + (1 + «) (g^(p 2 + q)^ + ^(Pl ~ qD) ■ (16) 

For k = —2 the vertex T ax reduces to the standard jWW-vertex. We can now express the 
matrix elements including the gauge link in the one gluon approximation in the following way 



(sdq; P-p\ W[oo,0,0r|0]</>i(0)|P,S) 



l-gl 



-le 



q&dq 



(2tt)< 



g sc ((l +p y)V SC (P-p-l) 



[ty + f + m q )u(P, 5)] i v(2P-2p- I) 
[l-v + iO] [I 2 + tO] [(I + p) 2 -m 2 + iO] ' 

(17) 



{adq; P - p, X\ W[oo, 0, T | 0] ^i(O) \P, S) 



l-gl 



dH 9ax{{p + l?) <{P-vA)v- { P-p-i), 



' qaq 1 (2tt) 4 ^3 ~ aK ~ "'v- w < 
[ g yp v . ( 2 p - 2p - I) + (1 + «) (t> CT (P - p + iy + ^ (P - p - 2l) a )] 
[l-v + iO] [I 2 + iO] [(I + p) 2 -m 2 q + iO] 

pr)\ 



(j* + / + m,) 75 7" 



where the subscript, 1 — gl denotes "one gluon exchange". In these expressions V denotes the 
propagator of the scalar and axial- vector diquark, 



V sc {P-p-l) 



-D^(P-p-l) 



[{P-p- I) 2 -m 2 + iO] ' 

/ _ (p- p -0.(p-p-0, 



(19) 



(20) 



[(P-p- I) 2 -m 2 + i0] 

The term ^ p p p j s a crucial difference of our approach compared to the calculation 

in Ref. [IT] , where the dependence on the proton and spectator momenta inside loop integral 
is absent. It is shown below that this leads to various complications when performing the 
loop-integral. 

In the similar fashion as for /i and h^ L , we extract the Boer-Mulders function by inserting 
Eqs. (fT7]) and (jT8j) (and the tree-matrix elements j5j) and ((11), i.e. the leading non-trivial 
perturbative contribution is the interference term between tree-graph and box-graph) into the 
quark-quark correlator (|3j) 

24p j T hj;{x,f T ) = ^ I dp- (Tr [$ unpo i(p, S)ia l+ ^} + Tr [$ unpo i(p, -S)ia i+ lb ] ) 

1 J p+=xP+ 

(21) 

where = e~ +y and e 0123 = +1. 
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Figure 2: Contribution of the gauge link in the one-gluon approximation. Left Panel: Box-graph. 
Right Panel: Box-graph hermitian conjugated. 



3 Boer-Mulders function for an axial-vector diquark 

We proceed with calculating the Boer-Mulders function in the axial- vector diquark sector. As 
described above, the interference term between tree- and box graph reads 



CqCdq 



1 



M 



dH 



8(2tt) 3 p| + m 2 P+ 
V m {P - P -l)(^2 4(P ~ P; A)^(P - p; A) J x 

[g a P v{2P-2p-l) + {\ + k) {v a (P-p + l) p + v p (P-p- 2l) a )] 

[l-v + iO] [I 2 - A 2 + iO] [(I + p) 2 -m 2 + iO] X 

(f + M) (V - R g ^pj (i> - m q ) 7 +f x (/ + j> + m q ) Cf + R g ^\ 75 



Tr 



The momentum of the quark, p is specified by 

p% + m 2 s - (1 -x)M 2 



P 



P 



2(1 -x)fH 



xP + , pr 



+ h.c. . 

(22) 

(23) 



3.1 Light cone integration 

Here we comment on the evaluation of the four-dimensional loop integral in Eq. [22j A conve- 
nient way to simplify the calculation is to sort the numerator in terms of loop momenta and 
consider each term separately. Since the numerator in Eq. (|22l) contains at most the loop 
momentum to the power of four we can write it in the following way 

4 

numerator = ^ N^ a .l ai ...l ai + N<®. (24) 
i=i 

The (real) coefficients (tensors) Na}... ai depend only on external momenta P (nucleon mo- 
mentum) and p (quark momentum) and can be computed in a straight-forward but tedious 
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calculation. We used the Mathematica-p&ck&ge TRACER [58] for this decomposition tf24|) 
The advantage of this procedure is that we are left with an arbitrary integral of the form 



(2tt) 4 [l-v + iO] [I 2 - A 2 + iO] [(/ + p - P) 2 -m 2 + iO] [(I + p) 2 -m 2 + iO] ' 

(25) 

and the light cone components of the loop momentum, l + and I , can be integrated out easily. 

We sketch the light cone integration. First, we specify the vector v to be a light cone 
vector v = [v~ = l,v + = 0,vt = 0] representing the Wilson line. Thus, the product I ■ v 
reduces to l + and doesn't contribute to the l~ integration. Next, we perform the integral over 
l~ via contour integration and encounter three poles in the Z~-plane from the last three terms 
in the denominator in Eq. (|25l) . The integral is non-vanishing when — xP + < l + < (1 — x)P + , 
otherwise all poles are located in the same complex l~ half-plane. For —xP + < l + < (1 — x)P + 
the third factor in the denominator in Eq. (|25l) has always a positive imaginary part while 
the forth factor has always a negative one. The imaginary part of the second factor becomes 
positive for l + < and negative for l + > 0. We close the contour of integration in the upper 
half-plane which excludes the forth factor in the denominator, and the second for Z + > 0. 
Thus, we obtain 

- r{l-*)P+ dl+ 1 

(26) 



T(*)e»i...a i _ _ • 

J ~ 3 



(27T) 



-xP+ 



2tt [l+ + i0] 
1 



21+ 2(1+ - (1 - x)P+)(P- -p~)+ [{l T + Pt) 2 + ' 



2(1+ - (1 - x)P+)\l 2 + X 2 



[2(1+ - (1 - x)P+)] [g ax ((l+P) 2 )9* ax (P 2 )(l ai -l a >)] 



I--P- 77- I ( r T+PT) 2 + m l 
' — P _P +2(I+-(l-x)P+) 



2(1+ - (1 - x)P+) (2(1+ + xP+)P- - [(It+Pt) 2 + m2j) + 2(1+ + xP+)\(T T +p T f 

[2i+}e(-i+)[ gax ((i + p) 2 )g: x (p 2 )(i ai -^)} 



+ mi 



21 + 



21+ (2(1+ + xP+)p- - [(Tt+Pt) 2 + m2]) + 2(1+ + xP+)[P T + A 2 ] 



We calculate the Z + -integral by adding the complex conjugated integral (stemming from the 
complex conjugated interference graph). Since l + + i0 is the only propagator remaining with an 
imaginary part, adding the complex conjugated integral results into a 5-function contribution 
via the relation 

' ' -2m6(l + ). (27) 



1+ + i0 1+ - iO 



We obtain 

j(i)ai— «i _|_ ^j(i)ai...ati\' 



(2f 



i r d 2 i T 

3 J Wf 



(l-x)P+ 



(-1) {g ax ((l+p) 2 )gt x {p 2 )[l a K..l^]) 



4_ ( r T+p - T) 2 

' - 2(l-x)P+ > l ° 



2P^ 



/| + A 2 



(l T +p T ) 2 + m 2 



dr 



cP+ 



1+ 

~xP+~ 



®(-l+)5(l+){g ax ((l + p) 2 )gl x (p 2 )rK..l^}) 



' 21+ 



2(1 - x)P+[P T + \ 2 } P T + \ 2 
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At this point we are forced to specify the form factor g ax since the second integral in 

Eq. (|28l) is potentially ill-defined. This happens when g(p 2 ) is a holomorphic function in 

p 2 (i.e. it contains no poles) and at least one of the Minkowski indices is light-like in the 

minus direction, e.g. a± = — , 0:2, •■•,0;? £ {+,_!_}. In such cases we end up with an integral 

of the form J dl + 5{l + )Q{— 1 + ), which is ill-defined. This implies that l + = and l~ = 00 

which signals the existence of a light cone divergence as was shown in Ref. [55]. While for a 

scalar diquark one doesn't encounter a Minkowski-index, otj = — , for twist-2 "T-odd" PDFs 

such as the Boer-Mulders function (so that there are no light cone divergences in this case), 

it was shown in Ref. [55] that for twist-3 "T-odd" PDFs light cone divergences exist for a 

u\ 

scalar diquark. However, calculating the coefficients Nai... ai in Eq. (I24|) for an axial-vector 
diquark, one of the Minkowski-indices can be a "minus" (i.e. otj = —). Thus, for an axial- 
vector diquark we encounter a light cone divergence already for twist-2 "T-odd" PDFs. Here it 
is worth mentioning that such divergences do not arise in a pQCD-quark-target model where 
the spectator state is a gluon |59|. 

From the standpoint of phenomenology one can regard these light cone divergences as 
model artifacts (for the axial- vector diquark model). It was shown in Ref. [55] how to regularize 
these light cone divergences by introducing non-light-like Wilson line. It was pointed out that 
one can also handle the light cone divergences by introducing phenomenological form factors 
with additional poles. Like the quark propagator in Eq. (|22l) they introduce additional factors 
of Z + in the numerator of the second term in Eq. ([28]) , We adopt this procedure to model the 
Boer-Mulders function and choose a form factor of the following form 



and find that for n > 3 (for n > 2 is already sufficient for k = — 2 ) enough powers of l + enter 
the numerator of the second term in Eq. ([281) to compensate the minus components l~. The 
second term then vanishes since J dl + l +n 5{l + )Q(—l + ) = for n > 1. In Eq. (|29p f(p 2 ) is 
then a function without poles to be fixed below, while A is an arbitrary mass scale to be fixed 
by phenomenology (i.e. fitting /1 to data). N ax is a normalization factor. 

After performing the light cone integrations we are then left with the remaining integral 
over the transverse loop momentum (for n = 3) 



9ax{p 2 ) = K 



ax 



(p 2 - m 2 )f(p 2 ) 
[p 2 - A 2 + i0] n 



(29) 



e^p J T h^ ax (x,p^) 



e g e dq A , 4 1 (l-z) 3 /(p 2 ) 



x 



8(27T)3 aX 3 [ p -| + ^2]3 m 4 



\4pt YJt? a p + 2 (It ■ Pt)^t b v + ^t c p + 2 ( /_ r • Pr)D P + E p \ + 
ep T [(fefAi + 2(T T ■ p T )P T Bi + %Ci + 2(f T • Pr)A + E x 
+e r T T T p s T \(Ai p % + 2l T ■ ptBi p ){Vt + p\) + E lv (t T + 2p l T )\ \ , ( 





(30) 



where the coefficients which are functions of x, M, m^p 2 -,, k and R g are given in the Appendix 
lAl rh\ replaces rh 2 by m 2 = xm 2 — x(l — x)M 2 + (l — x)A 2 . We point out that the vanishing of 
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the coefficient E p ensures that no IR- divergence appears in the transverse integral / fgflj) . This 
serves as an important check of the calculation. 



3.2 Transverse integral 

The transverse integral ([30]) can be calculated in a straightforward manner. We note that 
the transverse integral (|30j) is UV-divergent if we chose f(p 2 ) = 1. This can be seen from 
naive power counting of the integrand. The UV-divergence stems from the term Ij,, which in 
turn is a consequence of the fact that we took the full numerator of the axial- vector diquark 
propagator into account, in contrast to Ref. [41 J where no UV-divergences were reported. We 
regularize it by choosing f(p 2 ) to be a covariant Gaussian, 



M+pY 



-b\i+p\ 2 



(31) 



where b is interpreted as the the regulator of the high I integration. Due to the pole contribution 
l + = 0, the Gaussian has no effect on the light cone integration. Thus we can write the squared 
products of the momenta l+p and p - after performing the light cone integration - as follows, 



(i + pY 



P 



{It + Pt) 2 + xm 2 s - x(l - x)M 2 
1 — x 

p\ + xm 2 — x(l — x)M 2 
1 — x 



(32) 
(33) 



Now all that remains is to perform the It integration. After a shift of the integration variable 
from It — > It + PT it is convenient to use polar coordinates to calculate the integral, and 
perform the angular integration first. For this we choose a coordinate system in such a way 
that the x-axis is along pr, such that pr = |pr|(l, 0). The integration is performed with 
respect to that direction, i.e. It = |/t|(cos 4>, sin^>). Having fixed the coordinate system in 
such a way, the hanging index i can only be 2. We perform now the angular integration over 
4> by means of the formula 



n cos(nx)dx tt 
o l + acos(x) v 7 ! - a 2 



vT 



l 



a 2 < 1, n > 0. 



(34) 



We are left with the remaining one-dimensional integrals (y/z = I) 



hi 



,ax [ -2 \ 

(x,P T ) 




(1 - x) 3 e- b( Pr 



+1xm 2 a -2x(l-x)M 2 ) 



in, 



[?T + 



ml] 3 



[z(A p - 2At + Bt) + f T {A p -At- 2{B P - £,)) + C v - Q 



2 Alp 
z ^dr + Z 



2p 



-A 



1 p & 



+ 



Ei r _ 



P T ' 2 Pt 



2(D p - D t 



El 

-n 

P T 



^-? t ^ + d 1 + \e 1p 



(35) 
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Figure 3: Left Panel: The unpolarized u- and rf-quark distributions functions versus x compared to 
the low scale parameterization of the unpolarized u- and d- quark distributions [60]. Right Panel: The 
half moment of the Sivers functions and the unpolarized u and d distributions versus x compared to 
the low scale parameterization of the unpolarized u- and d- quark distributions (k = 1.0). 



where 6 = 6/(1 — x). These integrals can be expressed in terms of incomplete T-functions 

/>oo 

Y{n,x) = / e-H^dt. (36) 



The Boer-Mulders function for an axial-vector diquark then reads 



,ax [ -8 \ 

(x,P T ) 



e q e dq Ni (\-xfc b PT-' 2 ' 



48(2vr 



,4 - 1 " ax 



111 



Kt' ax [x,p^R g ,K,b,A,{M}) , 

(37) 



where the explicit form of H± is expressed in term of incomplete Gamma functions and can 
be found in Appendix O 

The Boer Mulders-function for a scalar diquark is much easier to calculate [38] due to its 
simpler Dirac-trace structure. The light-cone divergences we have encountered in the axial- 
vector diquark sector do not appear for the scalar sector. With our choice of the form factor 
(cf. Eqs. (|29j) and ([31]) ) the Boer-Mulders function for a scalar diquark reads 



h{ ,sc {x,p%) 



e Q e dq AT 4 (1 - xfM(xM + m q ) 



4(2vr) 



-N 

4 sc 



213 



-b(p2,+2xm 2 a -2x(l-x)M 2 ) 



p<(m brhD - m h(f T + mi)) + 



1 - b (j% + ml) 
2{f T + m\) 2 



(38) 



3.3 Sivers-function in the diquark spectator model 

Having obtained the results for the Boer-Mulders function hi, it is straight-forward to apply 
the procedure described above to calculate the Sivers-function The Sivers-function can 
be extracted from the following trace of the quark-quark correlator j3]) (see e.g. [30l 131"] ). 



M 



2S T e T ! p T ft T (x,p T ) = ^ dp' (Tr [ 7 +$(p; P, St)] - Tr [ 7 +$(p; P, -St)] ) 



■ (39) 
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Figure 4: Top Panels: The half-moment of the Boer Mulders (left) and Sivers (right) functions versus 
x compared to the unpolarized u- and d- quark distribution functions. Bottom Panels: The half 
moments of the Boer Mulders and Sivers functions, K = 1.0 lower-left, K = —0.333 lower-right, versus 
x, extractions from data were presented in Ref. [6T| 162. 63, 64, 65J for the Sivers function. 



It is well-known [48j that in the scalar diquark spectator sector the Boer-Mulders function 

and the Sivers function coincide, so the scalar Sivers function is given by the left-hand side of 

Eq. (|38l) . By contrast the different Dirac structure for the chiral even Sivers function and chiral 

odd Boer-Mulders function in the axial-vector diquark sector, Eq. (|21~T) and ([221) respectively, 

m 

lead to different coefficients in the decomposition Nai... ai in Eq. (1241) . So, whereas the form 
of the Sivers-function f^ ax is the same as the form of h^ ,ax given in Eq. (|37l) . the coefficients 
A p , B p , C p , D p , E p , Ai, B h Q, D h E h A lp , B lp , E lp differ. They are given for f^ ax explicitly 
in the Appendix [Bj 

4 Asymmetries in SIDIS 

4.1 Azimuthal cos(20)-asymmetry in unpolarized SIDIS 

Almost 30 years ago it was pointed out that both kinematic [66] and dynamical effects [67, 68J 
could give rise to an cos2</) azimuthal asymmetry going like p\/Q 2 (where Q is a hard scale) 
when transverse momentum scales are on the order of the intrinsic momentum scales of partons, 
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Figure 5: The first moment of the Boer-Mulders and Sivers functions versus x for k = 1.0. 
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Figure 6: Left Panel: The cos 20 asymmetry for 
kinematics. Right Panel: The cos 20 asymmetry 
kinematics. 



0.04 




7T + and 7r as a function of Pt at JLAB 12GeV 
for 7r + and 7r~ as a function Pt for HERMES 



Pt ~ pt- However, when transverse momentum is on the order of the hard scale, Pt ~ Q, 
these non-perturbative effects are expected to decrease relative to perturbative contributions 
[69, 70j. By contrast, taking into account the existence of "T-odd" TMDs and fragmentation 
functions it was pointed out by Boer and Mulders [8] that at leading twist a convolution of the 
Boer-Mulders and the Collins functions would give rise to non-trivial azimuthal asymmetries 
in unpolarized SIDIS. 

Having explored the flavor dependence of the hj~ we are now in a position to extend 
early phenomenological work on "T-odd" contributions to azimuthal asymmetries in SIDIS 
performed under the approximation of scalar diquark dominance|49|. In particular we consider 
the spin independent double "T-odd" cos 2(p asymmetry for tt + and tt~ production. 

The general form of the cross section for an unpolarized target reads [31 J 

da 2ira 2 \ ( _ 1 2 



dx dy dz d4>h dP? i xyQ 2 



1 - y + ) Fuu,t + (l - y) F UUtL (40) 



+ (2 - y) cos(^) F$Z+ + (l-y) cos(20 h ) F™ 2 ^ 
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Figure 7: Left Panel: The cos 2<f> asymmetry for ir + and ir~ as a function of x at JLAB 12GeV and 
HERMES kinematics. Right Panel: The cos 20 asymmetry for ir + and ir~ as a function z for JLAB 
kinematics. 



where the structure function F^ 2 ^ h is of most interest for the purpose of this paper. At 
leading twist it factorizes into a convolution of the Boer-Mulders and Collins fragmentation 
function pi [31] 



pcos 2<j> h _ p 

r uu — <- 



2h ■ k T h p T -k T p T . ± 



(41) 



MM h 

the convolution integral C is given by 
C[wfD] = xJ2^af d 2 p T d 2 k T 5^(p T -k T - P h± /z)w(p T ,k T ) f a (x,p 2 T ) D a (z,k 2 T ), (42) 

where summation runs over quarks and anti-quarks. p T , &t are the intrinsic transverse mo- 
menta of the active and fragmenting quarks respectively and Ph is the transverse momentum of 
the fragmenting hadron with respect to the photon momentum q. h is defined as Ph±/\Ph±\- 
We have fixed most of the model parameters such as masses and normalizations by com- 
paring the model result for the unpolarized "T-even" PDF f\ for u- and d- quarks (Eqs. ((H) 
and (fTUI) ) to the leading order (LO) low-scale {p? = 0.26GeV) data parameterization of Gliick, 
Reya, and Vogt [60]. Note that PDFs for u- and d- quarks are given by linear combinations 
of PDFs for an axial vector and scalar diquark, u = |/ sc + \ f ax and d = f ax [53l HTj. The 
best model approximation to the GRV data parameterization for u and d of [60J is shown in 
Fig. El and corresponds to a set of parameters m q = GeV, m s = 1.0 GeV, m ax = 1.3 GeV, 
A = 1.3 GeV, M = 0.94 GeV - fixed, and R g = 5/4. For "T-odd" PDFs such a procedure 
for fixing the model parameters is not sufficient since it doesn't determine the sign and the 
strength of the final state interactions. In our case the final state interactions are described 
effectively in the one gluon exchange approximation by the product ed q e q , the charges of the 
diquark and quark, respectively. We need to fix the value of this product. For that reason we 
calculated the Sivers function for u- and d- quarks in the diquark model and compared our 
results in Fig. [3] for the "one-half" moment, 

/i L T (9 ' 1/2) W= /dV§Ur (,) (*,P§0. (43) 
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as well as the first moment, 



ft^ (*) = / d 2 PT ^ f#<> (s, &) , (44) 



with the existing data parameterizations where q represents the quark flavor (see Refs. [61, 
1621 1631 1641 [65]). In such a way we are able to fix edqe q /47r = Cpa s = 0.267 with color 
factor, Cf = 4/3. We display the "one-half" and first moments for u- and d- quark Sivers 
functions f^p and Boer-Mulders functions hf^ (where q = u,d) along with the unpolarized 
u- and d- quark pdfs in Figs. [4] and The "one-half" and first moments of the u- and 
d- quark Sivers functions are negative and positive respectively while the u- and d- quark 
Boer Mulders functions are both negative over the full range in Bjorken-x. These results are 
in agreement with the large N c predictions [42J, Bag Model results reported in [40J, impact 
parameter distortion picture of Burkardt [45] and recent studies of nucleon transverse spin 
structure in lattice QCD [H]. H Also, we explored the relative dependence of the d-quark 
to u-quark Sivers function, see Fig. [H For example, choosing a value of k = —0.333 as was 
determined in Ref. [57] we find the d quark Sivers is smaller than the u- quark. Choosing k = 1 
we find reasonable agreement with extractions reported in [61]. It is worth noting (see Fig. [4] 
) that the resulting u-quark Sivers function and Boer-Mulders function are nearly equal, even 
with the inclusion of the axial vector spectator diquark. An exact equality was first noted in 
the simpler scalar di-quark dominance approximation in [48] . 

Our model input for the Collins functions is based on very recent work in |71] where the 
Collins function was calculated in the spectator framework. Therein it was assumed that 
jj±.(dts fav) ^ _jj±(fav) ^ e pj on sec t r, thereby satisfying the Schaefer-Teryaev sum rule 



2} locally. We use those results along with the results of this paper for to estimate the 
azimuthal asymmetry A^y 2 ^ (cf. Eq. ([41~]) ). where 

a cos 2$ _ Jjg cos2 ^ d(T 

Auu = fd$da (45) 

and d& is short-hand notation for the phase space integration. In Fig. [6] we display the 
A c ° b v 2(t> (P T ) in the range of future JLAB kinematics [73] (0.08 < x < 0.7, 0.2 < y < 0.9, 
0.3 < z < 0.8, Q 2 > lGeV/c, and 1 < E v < 9 GeV) and HERMES kinematics [1] (0.23 < 
x < 0.4, 0.1 < y < 0.85, 0.2 < z < 0.7, with Q 2 > 1 GeV/c, and 4.5 < E w < 13.5 GeV). In 
Fig. [7] we display the x and z dependence in the range 0.5 < Pt < 1.5 GeV/c. It should be 
noted that this asymmetry was measured at HERA by ZEUS, but at very low x and very high 
Q 2 [35J where other QCD effects dominate. It was also measured at CERN by EMC [74], but 
with low precision. Those data were approximated by Barone, Lu and Ma [75] in a w-quark 
dominating model for h^, with a Gaussian, algebraic form and a Gaussian ansatz for the 
Collins function. Our dynamical approach leads to different predictions for the forthcoming 
JLab data. 



lr rhis is in agreement with the value a a used in [71] to explore the "T-odd" fragmentation functions. It is 
worth noting that the running coupling extrapolated to the "low-scale" fi 2 in [60j is different than the coupling 
that characterizes the FSIs in the one gluon exchange approximation. 

2 It is interesting to note the approximate agreement for the flavor dependence of h^t among such models 
probably arises because our input quark and di-quark wave functions share the same SU(4) flavor-spin de- 
pendence as the bag and other spectator models. Additionally SU(4) symmetric baryon wave functions are 
compatible with large- iV c counting rules. 
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Figure 8: Left Panel: The half-moment of xhf^ 2 ^ versus x compared to the unpolarized u- and d- 
quark distribution functions. Right Panel: The sin 2(f) asymmetry for n + and tt~ as a function of x at 
JLAB 12GeV kinematics. 



4.2 Single spin asymmetry A^f 2 ^ in SIDIS 

Since we have calculated the chiral-odd but "T-even" parton distribution h^ L (cf. Eqs. (fl~3j) . 
(|T2l) ) we use this result together with the result of Ref. [71 J for the Collins function to give 
a prediction for the sin(2</>) moment of the single spin asymmetry Ajjl for a longitudinally 
polarized target. In particular we are able to take into account the flavor dependence of the 
asymmetry. We adopt the similar procedure for the azimuthal cos(2(/>)-asymmetry for treating 
the leading twist observable A s ™^ 2 ^ . 

A decomposition into structure functions of the cross section of semi-inclusive DIS for a 
longitudinally polarized target reads (see e.g. [3T] ) 



dauL 



dx dy dz d(ph dP, 



h± 



2vra 2 
xyQ< 



Sii 



(1 - y) sin(2«^) F*%W + (2 - y)7r^ S in( 



UL 



(46) 



where S» is the projection of the spin vector on the direction of the virtual photon. In a 

partonic picture the structure function F^ 2 ^ jg a leading twist object (while is sub- 

leading) and given by a convolution of the TMD h\ L and the Collins function (cf. |31| ) 



F, 



UL 



c 



2h ■ kx h ■ p T — kx • Pt 

MM h 



htr.H 



(47) 



where the explicit form of the convolution is given in Eq. (|42l) . 

We insert our result for h{ L (Eqs. QUI), (fl^l ) and the result of Ref. [71] into Eq. gZJ 
to compute the single spin asymmetry. This is the first calculation of this observable in the 
spectator framework, whereas the part of F^f^ described by higher twist "T-odd" PDFs has 
been analyzed in the diquark model in Refs. [2U EHJ [28] . Similar phenomenology for F^ 2 ^ 



and F^^ has been performed in Refs. [761E7] using the framework of the chiral quark soliton 



model. 
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We display the results for the single spin asymmetry A^£ in Fig.[8]using the kinematics 
of the upcoming JLab 12 GeV upgrade. We note that the tt~ asymmetry is large and positive 
due to the model assumption H^ dls f av ^ ~ _j{^-(f av \ This asymmetry has been measured 
at HERMES for longitudinally polarized protons [TS] and deuterons [21]. The data show that 
for the proton target and HERMES 27.5 GeV kinematics both ir + and ir~ asymmetries are 
consistent with down to a sensitivity of about 0.01. That is to say, these asymmetries could 
be non-zero, but with magnitudes less than 0.01 or 0.02. These results are considerably smaller 
than our predictions for the JLab upgrade. For the deuteron target the results are consistent 
with for 7r + and tt~ . There is one ir° point at x ~ 0.2 that could be positive at about 0.03. 
This SIDIS data for polarized deuterons could reflect the near cancellation of u- and ci-quark 
h^ L functions and/or the large unfavored Collins function contributions. There is also CLAS 
preliminary data [78] at 5.7 GeV that shows slightly negative asymmetries for ir + and ir~ and 
leads to the extraction of a negative xh^ . This suggests that the unfavored Collins function 
(for d — ► 7r + ) is not contributing much here, unlike the inference from the HERMES data. 
Data from the upgrade should help resolve these phenomenological questions. 



5 Conclusions 

In this paper we performed calculations of transverse momentum dependent parton distri- 
butions, including the Boer-Mulders function hf, the Sivers function and h\ L in the 
diquark spectator model taking into account both axial-vector and scalar contributions. The 
calculation of these functions in both sectors allowed us to explore their flavor dependence, 
i.e. to compute their u-quark and a d-quark contributions. For "T-even" distributions like 
h± L , a non-trivial contribution could already have been obtained from a tree-level diagram. 
By contrast, final state interactions or, equivalently, contributions from the gauge link had 
to be taken into account for the "T-odd" Boer-Mulders and Sivers functions, requiring the 
calculation of a loop (box) diagram. It was found that the loop integrals for f^f and show 
light cone divergences and UV divergences in the axial-vector diquark sector, in contrast to 
the scalar diquark sector. We regularized these divergences by choosing specific types of the 
phenomenologically motivated nucleon-diquark-quark form factors. By comparing the model 
expression for the unpolarized parton distribution, f\, with the low scale parameterizations of 
that function obtained from data, it is possible to fix most of the parameters of the model, 
masses, normalization and R g , the ratio of axial diquark couplings to the nucleon. In order 
to fix the sign and size of the final state interactions specific for "T-odd" distributions, we 
calculated the Sivers function and compared the result to parameterizations of SIDIS data 
already determined in a global fit for fy^. In such a way the remaining parameters could be 
fixed, and predictions were presented for the flavor dependence of the Boer-Mulders function 
hj^ for a u-quark and a d-quark. We find that the fop-half- and first-moments of this function 
are negative for both flavors. This result is in contrast to that in [41] . Our sign result is in 
agreement with other approaches that predict negative hf and hf^ d \ 

We also note that the resulting u-quark Sivers function and Boer-Mulders function are 
nearly equal, even with the inclusion of the axial vector spectator diquark. This near equality 
h^ ~ f^p was obtained from models without axial di-quarks, hinting at some more general 
mechanism that preserves the relation. 

We used our result for h^ as one ingredient in the factorized formula for the azimuthal 
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asymmetry A^^ in unpolarized semi-inclusive lepto-production of positively and nega- 
tively charged pions. We also used our h\ L as an ingredient in the single-spin asymmetry 

Ay 1 ^ 2 ^ for a longitudinally polarized target in semi-inclusive DIS. Another key ingredient 
for determining such asymmetries is the Collins fragmentation function . For this function 
we used the most current expressions that were obtained in a similar spectator model. 

We provide estimates for A^^ 2 ^ and A^ 1 ^ 2 ^ . The latter has already been measured at 
HERMES and preliminarily by CLAS. There are important differences in the kinematic regions 
explored, but there remain discrepancies that may be resolved in the future at Jefferson Lab, 
for which our model gives striking predictions of relatively large asymmetries. The non-trivial 
7r~ asymmetry is driven in large part by the model assumption ff^ dls f av ^ ~ _jjHf av ) - m 

the pion sector. We note that our result for A^ 1 ^ 2 ^ is the first phenomenological treatment 
in the spectator framework of this observable. 

The former unpolarized asymmetry, A^^ 2 ^\ was measured at HERA, but for very small 
x and extremely high Q 2 . Again, this will be measured in the future at JLab. We predict that 
those results should correspond to the opposite sign asymmetries for opposite charged pions. 

In summary, a refined diquark spectator model, including axial vector di-quarks leads to 
both u- and d-quark "T-odd" TMDs and provides the ingredients for predicting a range of 
asymmetries for future experiments. The approach we have been taking is to use and refine a 
model for the soft regime that makes sense in QCD and can be applied broadly to a range of 
measurable phenomena. We have fixed the parameters in the model to approach the inferred 
structure of the lowest order asymmetries. Combined with the recent determination of the 
fragmentation functions, we have predicted new SIDIS results. The spirit of this work is to 
understand the dynamics of processes like SIDIS by refining a robust and flexible model for 
the "T-odd" functions that compares with existing data. While a global fit to all the data 
eventually can be performed, the underlying mechanism is likely to be revealed by honing in 
on more sophisticated and inclusive models, as we have done here. 
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A Appendix: Axial Diquark Contribution to Boer Mulders 
Function 

The coefficients appearing in the calculation of the Boer Mulders function for the axial vector 
diquark (see text) read 



A p = (1-x) (2Mm q + (2-R g )M 2 + R g m 2 )(3 + 2K) + R g m 2 s (7 + 4 K ) 



(48) 
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B p = (1-x) (2Mm q + (2-R g )M 2 + R g m 2 ){3 + 2K) + R g m 2 s (n + 6K) 



(49) 



C p = 2m 2 s (l-x) 2Mm q (-l + 2x(2 + K))+R g (m; + 6&(2 + K)+m 2 q (3 + 2K)) 
+M 2 (-2 + Ax(2 + k) + R g {\ - 2x(2 -x){2 + k))) 



(50) 



Dp = 2m 2 s (l-x) 2Mm q (-l + 2x(2 + k)) + R g (m 2 s + 2p 2 T {2 + k) + m 2 (3 + 2k)) 
+M 2 (-2 + 4z(2 + k) + R g (l - 2x(2 - x){2 + «))) 



(51) 



E p = 0, 



(52) 



A, = 



R g {3 + 2 K ) 
2M 



R g m q (m 2 s - v T ) + M 2 m q {2 - R g ){\ - x) z + x(l - xfRgM 3 



+M ((2 - R g )xm 2 -f T {2- xR g )) 



(53) 



B, 



Rg_ 

2M 



+m q R g -f T + m; 



(3 + 2k) [M 2 m q (2 - R g )(l - xf + x(l - x) 2 R g M 3 
2 7 + 4k\ 



3 + 2k J ^ 



+M[~f T {2-xR g )(i + 2K)+m 2 s {A{2 + K)+x(&-7R g +AK{l-R g )))] , (54) 



Cl = -7TTT 



1 

2M 



(1 - .x) 3 (3 + 2K)M A m q (-4 + (4 - i? s )(l - x)R g ) 



-(1 - x) 3 a;(3 + 2K)M b R g {2 - (1 - x)i? s ) 

+ro,i^ (8m^(2 + k) + (w 2 - f T ){m 2 s +f T ){2> + 2k)) 

+MR g {xm i s {2 -R g )(3 + 2K)-p^{2- xR g )(3 + 2k) 

-2m 2 s f T (-17 + x + 8xR g - 8k + AxR 9 k) + 2m 2 {\ - x)(m 2 s + (3 + 2k)p%^) 
+2(1 - x)M 2 m q (f T (3 + 2k) (2 + (1 - x)R g (2 - R g j) 

+m 2 s (A(2 + k) - x(6 + 4k))) + 2(1 - x)M 3 ( - m 2 q (2 - R g )(l - x) 2 (3 + 2k) 
+f T (2-R g + x(l - x)R 2 g ) (3 + 2k) + xm 2 s (2 + R g (6 - 7x + 4/c(l - a;)) )) 



(55) 



M 



m q R 2 (m z s +p z T ) + (1 - x) l M z m q R g {2 - R g ) 
+M 3 (1 - x) 2 (xR 2 g + 4(2 + k) - 2R g {\ + x){2 + k)) 
+MR g (m 2 s {2n - x(2 - R g + 2k)) + p%{x(4 + R g + 2k) - 2(3 + 



(56) 
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m q R 2 g {f T + m 2 s f + M 4 m q (l - x) 3 ( - 4 + (1 - a;) (4 - i? g )i? s 

+8x(2 + k)) - x(l - x) 3 M 5 R g (2 - (1 - x)R g - Ax{2 + k)) 

-2M 2 m q {\ - x) (m 2 (4 + (1 - x){2 - R g )R g - 2x(3 + 2k)) 

+f T (-2 - (1 - x){2 - R g )R g ) + Ax(2 + «)) ) 

+MR g ( - x(2 - R g )m 4 s + 2m 2 s f T (2k - 1 - z(l - i? g + 2k)) 

+2m 2 q (l - x) (m 2 Jl + 2k) - f T {3 + 2k)) + f T (x(8 + R g + 4k) - 2(5 + 2k)) ) 

-2(1 - x)M 3 (^ - m 2 q {l - xf{2 - R g ){3 + 2k) 

+xm 2 s (6 - 8a; + R g {-2 + R g + 3x- xR g ) - 2xk(2 - R g )) 

+f T ( - a;(l - x)R 2 g + R 9 (5 + 2k- Ax{2 + k)) + 2(-5 - 2k + 4a;(2 + k)) 



Ai p ^Bi p = -AR g m 2 s {l-x){2 + K), 



Ei p = A{l-x)rn 2 s (2Mm q {l-x){2 + K) + M 2 {l-x){2-R g {l-x)){2 + K) 
+R g (Km 2 s - f T {2 + k) 



B Axial Diquark Contribution to Sivers Function 

For the Sivers function the corresponding coefficients read 



7 + 4k\ 



^p iv ' = -(1-.t)(3 + 2k) (2Mm q + (2-R g )M 2 + R g m 2 q + R g m 2 s j^-\ 



D 



Siv. 



-(1 - x){3 + 2k) ( 2Mm q + (2 - R g )M 2 + R g m 2 q + R g m 2 s 



: 11 + 6k\ 
3 + 2k ) ' 



a 



Siv. 



-2{\-x)m\ 



2Mm q {-\ + 2x(2 + k)) + R g (m 2 s + Qp 2 T {2 + k) + m 2 q (3 + 2k)) 



+M 2 (-2 + 4a;(2 + k) + R g (1 - 2a;(2 - x){2 + k))) 



D 



Siv. 



-2(1 -x)m 2 s 



2Mm q {-\ + 2x{2 + k)) + R g (m 2 s + 2f T {2 + k) + m 2 q (3 + 2k)) 



+M 2 (-2 + 4x(2 + K )+R g (l- 2x(2 - x){2 + «))) 



£f v - = 0, 
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3_+2k 
2M Hg 



R g m q (m 2 s - p 2 T ) + M 2 m q (2 - R g )(l - x) 2 + x(l - x) 2 R g M i 



+M (xm 2 s (2 - R g ) - f T (2 - xR g )) 



m 



Siv. 



Rg_ 

2M 



(3 + 2n)(M 2 m q {2 - R g ){\ - xf + M 3 R g x(l - xf 



+m q R g (-f T + m 2 s ^j^j - Mf T {2 - xR g )) 
+Mm 2 s (4(2 + k) + x (6 - 7R g + 4/s(l - Rg))) 



1 

2M 



(1 - .T)' j M 4 m g (4 - (4 - Rg)R g {l - x)){3 + 2k) 



-x(l - x) 3 M 5 R g (2 - (1 - x)R g )(3 + 2k) 

+m q R 2 g(8my T (2 + k) + (m 2 - f T )(m 2 s + f T )(3 + 2k)) 

+MR g (xm 4 s {2 - R g )(3 + 2k) - p%{2 - xR g )(3 + 2k) 

+2m 2 s p 2 T (17 - x - 8xR g + 8k - AxR g n) + 2m 2 (l - x)(m 2 s + p|(3 + 2k))) 

-2M 2 m q {\ -x)(- f T (2 + (2 - Rg)Rg(l - x))(3 + 2k) 

-m 2 s (A(2 + K) -x(6 + 4k))) + 2M 3 (1 - x)( - m 2 q {2 - R g )(l - x) 2 {3 + 2k) 

+f T (2-R g + x(l - x)R 2 g ) (3 + 2k) + xm 2 s (2 + R g (6 - 7x + 4k(1 - a;))) ) 



' M 



- m q R 2 (mj + Vt) + m q M 2 {l - x) 2 R g (2 - R g ) 
+ (1 - x) 2 M 3 (xR 2 g + 4(2 + k) - 2R g {l + x)(2 + «)) 
+MR g (m 2 s (x(-2 + R g - 2k) + 2k) + $,(-2(3 + k) + x{A + R g + 2k))) 



' M 



- m q (m 2 s + f T ) R\ 



+M 4 m q (l - x) 3 (-4 + (4 - Rg)R g {\ - x) + 8x{2 + «)) 
-M 5 i? 9 a;(l - x) 3 (2 - (1 - x)R a - Ax{2 + k)) 

+2M 2 m q {\ - x) (f T (2 + (2 - R g )R g {l - x) - Ax{2 + «)) 
+m 2 s (2x + Axk - (2 - Rg)R g (l - a;))) 

+M# r ( - am* (2 - J? g ) + 2m^ (2k - 1 - ir(l - i? s + 2k)) 
+2m^(l - x) (m 2 s (l + 2k) - p%(3 + 2k)) 
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+^(-2(5 + 2k) + x(8 + R g + 4k))) 

-2M 3 (1 -x)[-m\{2- R g ){\ - x) 2 (3 + 2k) 

+xm 2 s (2 -Ax- R g (2 - R g - 3x + xR g ) - 2(2 - R g )xn) 

+p% (-x(l - x)R 2 g + R g (5 + 2k- Ax(2 + k)) - 2(5 + 2k - Ax{2 - 



zlSiv. 
A lk 



B. 



Siv. 
Ik 



E, 



Siv. 
Ik 



0. 



(69) 



(70) 



C ^-functions 

The 7£-functions appearing in the text are defined in the following. The 7£-function for the 
unpolarized PDF f\ for an axial vector diquark reads 



U'r [x,4;R g ,{M}) 



= M 



(p z T Y + 2(1 - x(l - x))m 2 p z T + x 2 m\ + 6x(l - x) 2 m q Mm 



+ (1 - x) 2 M 2 {f T + 2x 2 m 2 + (1 - x) 2 m 2 q ) + m 2 q (l - x) 2 (f T + 2m 2 s ) 

+R g M 2 ( - (f T ) 2 - x 2 m\ + (1 - (4 - x)x)m 2 s f T 

-(1 - x) 2 m 2 q (f T - m 2 s ) -M 2 (l- x) 2 (f T - x 2 m 2 s + (1 - xfrn 2 ^ 

+m q M (x{f T + m 2 ) 2 + 2x(l - x) 2 M 2 (f T - m 2 s + x(l - x) 4 M 2 )) 



+Rl 



+ (m s - (1 - x)M) 2 ) (f T + (xM - m q ) 2 ) 



(f T + (m s + (1 - x)Mf) 



The corresponding 7£-function for hf^ x reads 



n 



-L.ax 



1L 

= -m q R 2 g (m 2 s + f T ) 2 + 2m q M 2 {\ - x)[f T \2 + R g (l - x)(2 - R g ) 
+m 2 s [2(2-x)-Rg(l-x)(2-R g )}) 

-M A m q (\ - xf [4 - (4 - Rg)R g (l - x)] - x(l - x) 3 M 5 R g {2 - (1 - x)R g ) 
+MR g {m 2 s +f T ) [2(1 - x)m 2 q - xm 2 s (2 - R g ) - f T {2 - xR g )] 

+2(1 - x)M 3 ( - m\{\ - x) 2 {2 - R g ) + xm 2 s (2 + R g (2 - (1 - x)R g - 3a;)) 
+f T (2-R g + x{l-x)R 2 g }. 



(71) 



(72) 



(73) 



(74) 



The 7£-function for the "T-odd" functions hi and f^f is more complicated and contains in- 
complete Gamma functions (see text) with a = bfh\, c = b(p^ + m\), d = bp^, 



K 



f' ax (x,f T ;R g ,K,b,A,{M}) 



(75) 
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n^ ax (x,f T ;R g ,n~b,k,{M}) 



= (A p - 2A t + B{) 



■Se-( 1 + |,r ( o,«, + ^ 



+ \A„ - A, - 2(B P - B,)\ + C p - Q) 



Alp 

2p| 



(1 + 2a + y )e a (T(0, a) - T(0, c)) - — - e~ d - 



2c 2 



Ai L _£ i E± 
-6e Q (l + |)(T(0, a) T(0, c )) + ^-±£ - ^ C - d (2c - a(l - c)) 



f ( -2(D p -A)-$ 



y e a (r(o, a )-r(o, c )) + fr 



+ - 



-6e«(l + ^)r(0, C ) + A (2c - a (l- c ))c- d 



1 



(76) 
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